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Abstract 

We discuss certain issues related to the limitations of density parameters for a 
"radiation" -like contribution to the Friedmann equation using kinematical or ge- 
ometrical measurements. We analyse the observational constraint of a negative 
(1 + z) 4 -type contribution in cosmological models. We argue that it is not possible 
to determine the energy densities of individual components of matter scaling like 
radiation from astronomical observations. We find three different interpretations of 
the presence of the radiation term: 1) the FRW universe filled with a massless scalar 
field in a quantum regime (the Casimir effect), 2) the Friedmann- Robertson- Walker 
(FRW) model in the Randall Sundrum scenario with dark radiation, 3) the cosmo- 
logical model with global rotation. From supernovae type la (SNIa) data, Fanaroff- 
Riley type lib (FRIIb) radio galaxy (RG) data, baryon oscillation peak and cosmic 
microwave background radiation (CMBR) observations we obtain bounds for the 
negative radiation-like term. A small negative contribution of dark radiation can 
reconcile the tension in nucleosynthesis and remove also the disagreement between 
Hq values obtained from both SNIa and the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite data. 



1 Introduction 



In constraining the negative energy contribution to the Friedmann equation, a 
variety of astronomical observations were used, such as SNIa data [1,2], FRIIb 
RG data [3], baryon oscillation peak and CMBR observations. Although we 
can obtain the bounds on the sum of density parameters for all fluids scaling 
like radiation, it is impossible to determine the contribution coming from each 
component. In particular it is impossible to determine the negative energy 
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density arising from the Casimir effect or dark radiation and positive energy 
density arising from radiation because they scale with respect to redshift in 
the same way in the Friedmann equation. We note that the CMBR and big- 
bang nucleosynthesis (BBN) offer stringent conditions on this term, which 
can be regarded as an established upper limit on any individual components 
of negative energy density, and therefore on the Casimir effect, global rotation, 
discreteness of space following loop quantum gravity or brane dark radiation. 
We argue that even if the precise value of the density parameter for fictitious 
fluid is known from observations, it is still not possible to determine the nature 
(or, say, the origin) of the global rotation of the Universe using measurements 
basing on geometrical optics only, i.e. null geodesies ds 2 = 0. For example, if 
we consider the origin of this negative radiation-like term from global rotation 
effects, then it is not possible to come up with different forms of contribution 
leading to the same expression of — O r (1 + z) 4 type. There are, in principle, 
three different interpretations of the presence of the negative radiation term 
in the H 2 (z) relation, as follows. 

The Casimir effect is a simple observational consequence of the existence of 
quantum fluctuation [4]. The Casimir force between conducting plates leads 
to a repulsive force, like the positive cosmological constant. Moreover, dif- 
ferent laboratory experiments were designed to measure the Casimir affect 
with increased precision and thus strengthen the constraints on corrections 
to Newtonian gravitational law [5]. For a survey of recently obtained results 
in Casimir-energy studies see [6]. For our purpose, it is important that the 
effect of Casimir energy, which scales like radiation, can contribute into the 
H 2 (z) relation — crucial for any kinematical test. It is also interesting that 
the same type of contribution to the effective energy density can be produced 
by loop quantum theory effects in semi-classical quantum cosmology [7,8]. 
These effects give rise to an evolutionary scenario in which the initial singu- 
larity is replaced by a bounce. It is worth mentioning that the Casimir-type 
contributions arising from tachyon condensation are possible [9]. 

In the brane world scenario, our universe is a submanifold which is embed- 
ded in a higher- dimensional spacetime, called 'bulk space'. In the Randall 
and Sundrum scenario [10], the Einstein equations restricted to the brane 
reduce to a generalization of the FRW equation with additional terms. One 
of these terms, called dark radiation, is of considerable interest because it 
scales like radiation [11]. This term arises from the non- vanishing electric 
part of the five dimensional Weyl tensor. Dark radiation strongly effects both 
BBN and CMBR. Ichiki et al. [12] gave limits on the possible contribution as 
-1.23 < pdr/p^ < 0.11 from BBN and -0.41 < p dr /p 7 < 0.105 at the 95% 
confidence level from CMBR measurements. Let us note that a small nega- 
tive contribution of dark radiation can also reconcile the tension between the 
observed 4 He and D abundances [13]. 



2 



Let us consider Newtonian cosmology following Senovilla et al. [14]. Then we 
can define, following the authors, a homogeneous Newtonian cosmology with p 
and p having no spatial dependence, i.e. p = p(t) and p = pit), and we require 
that the velocity vector fields depends linearly on the spatial coordinates. 
Then the equation, which represents shear-free Newtonian cosmologies with 
expansion and rotation, well known as the Heckmann-Schiicking model [15] 
takes the following form 

2 p(t ) 2u 2 

where u is the rotation scalar and C is an integration constant. We inter- 
prete it in terms of curvature constant although in the Newtonian spacetime 
the curvature is zero. For our purpose, it is important that the effect of rota- 
tion produces a negative term in the Newtonian analogue of the Friedmann 
equation. 

In the Newtonian cosmology in contrast to general relativity rotation can 
appear in homogenous and isotropic space. In general relativity the effect of 
rotation are strictly related to non-vanishing shear. The homogeneous universe 
with non- vanishing shear may expand and rotate relative to local gyroscopes. 
The relation between the rotation of the universe and the origin of the rotation 
of galaxies was also investigated [16,17,18,19]. Additionally, the role of rotation 
of objects in the Universe, their significance and astronomical measurements 
was recently addressed by [20,21]). 



2 Observational constraints on the FRW model parameter with a 
negative radiation-like term 



Usually the fundamental test of a cosmological model is based on the lumi- 
nosity distance as a function of redshift [22]. For the distant SNIa, one 
can directly observe their apparent magnitude m and redshift z. Because the 
absolute magnitude M. of the supernovae is related to its absolute luminosity 
L, the relation between distance modulus p, the luminosity distance cIl, the 
observed magnitude m and the absolute magnitude M has the following form 

/i = m-M = 51og 10 d L + 25 = 51og 10 £> L + Al (2) 

where D L = H cIl and A4 = —5\og 10 H + 25. The luminosity distance of a 
supernova can be obtain as the function of redshift: 
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where 

(|P) =fi m ,o(l + ^) 3 + ^,o(l + ^) 4 + ^r,o(l + ^) 4 + fidr,o(l + ^) 4 + fiA,0, (4) 

^fc,o = — an d J~{ x ) = (sinh(x), x, sin(x)) for k < 0, k = 0, k > 0, respec- 
tively. We assumed fi r o = fi 7 o+^o = 2.48/i~ 2 x 10~ 5 + 1.7/i- 2 x KT 5 ~ 0.0001 
[11]. 

Daly and Djorgovski [23] (see also [24,25]) proposed the inclusion of radio 
galaxies in the analysis. To accomplish this, it is useful to use the coordinate 
distance defined as: 

yv ; c(l + z) K ' 

Daly and Djorgovski [3] have compiled a sample comprising the data on y(z) 
for 157 SNIa in the Riess et al. [1] Gold dataset and 20 FRIIb radio galaxies. 
In our data sets we also include 115 SNIa compiled by Astier et al. [2]. 

The coordinate distance y(z) does not depend on the value of H from defini- 
tion. However, when we compute the coordinate distance from the luminosity 
distance (or the distance modulus /z) the knowledge of the value of Ho is re- 
quired. For each sample we choose the values of H which were used in original 
papers. We used the distance modulus presented in Ref. [1,2] for the calcula- 
tion of the coordinate distance. For each sample we choose the values of H 
apropriate to the data sets. For Riess et al.'s Gold sample we fit the value 
of h = 0.646 as the best fitted value and this value is used for calculation 
of coordinate distance for SNIa belonging to this sample. In turn the value 
h = 0.70 was assumed in the calculations of the coordinate distance for SNIa 
belonging to Astier et al.'s sample, because the distance moduli /i presented 
in Ref. [2, Tab. 8] was calculated with such an arbitrary value of h = 0.70. 
The error of the coordinate distance can be computed as 

where <Tj(/Zj) denotes the statistical error of distance modulus determination 
(note that for Astier et al.'s sample the intrinsic dispersion was also included) 
and a(H ) = 0.8 km/s Mpc denotes error in H measurements. 

Recently Eisenstein et al. have analysed the baryon oscillation peaks (BOP) 
detected in the Sloan Digital Sky Survey (SDSS) Luminosity Red Galaxies 
[27]. They found that value of A 




A = ~~ ~t , ■ ■HVlftfc.ol/. (7) 
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(where E(z) = H(z)/H and z x = 0.35) is equal A = 0.469 ± 0.017. The 
quoted uncertainty corresponds to one standard deviation, where a Gaussian 
probability distribution has been assumed. 

It is also possible include in our analysis the so called the (CMBR) "shift 
parameter" 



where R = 1.716 ± 0.062 [28]. 

In our combined analysis, we can obtain a best fit model by minimizing the 
pseudo-x 2 merit function [29]: 



where A 1 and R mod denotes value of A and R obtained for particular set of 
the model parameter. For Astier et al. SNIa [2] sample additional error in z 
measurements were taken into account. Here aiiyji) denotes the statistical error 
(including error in z measurements) of the coordinate distance determination. 

In order to constrain the cosmological parameters, one can minimize the fol- 
lowing likelihood function £ oc exp(— x 2 /2). However, to constrain a particular 
model parameter, the likelihood function marginalized over the remaining pa- 
rameters of the model should be considered [29]. Our results are presented 
in Table 1, Table 2 and Fig. 1. Table 1 refers to the minimum x 2 method, 
whereas Table 2 shows the results from the marginalized likelihood analysis. 

We obtain as the best fit a flat universe with f2 mi0 = 0.3, f2 dr ,o = and 
^a,o — 0.7. For the dark radiation term, we obtain the stringent bound f2 dri0 > 
-0.00035 at the 95% confidence level (-0.00026 at the 68.3% confidence level). 

Please note that if Q Tj0 + fi dr < 0, then we obtain a bouncing scenario 
[30,31,26] instead of a big bang. For fi mj0 = 0.3, fi dr>0 = -0.00035 and h = 0.65 
bounces (H 2 = 0) appear for z ~ 1200. In this case, the BBN epoch never 
occurs and all BBN predictions would be lost. 

To obtain stronger constraints on the model parameters, it is useful to use 
the CMBR observations. The hotter and colder spots in the CMBR can be 
interpreted as acoustic oscillations in the primeval plasma during the last 
scattering. It is very interesting that the locations of these peaks are very 
sensitive to variations in the model parameters. Therefore, they can be used 




X — XSN+RG + XSDSS + XCMBR _ 
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as another way of constraining the parameters of cosmological models. The 
acoustic scale £a which gives the locations of the peaks is defined as 



7T 



fZde 

Jo 

J Zdec 6 



dz' 
H(z') 
dz' 
H(z') 



(9) 



where, for the flat model, equation (4) reduces to 

h(z) = H oy Jn mfi (i + zf + a,o(i + ^) 4 + Mi + zf + fi A)0 , (10) 

where c s is the speed of sound in plasma. Knowing the acoustic scale we 
can determine the location of the m-th peak l m ~ Ia[vh — <f) m ) where <p m is 
the phase shift caused by the plasma driving effect. The CMBR temperature 
angular power spectrum provides the locations of the first two peaks l\ = 
220.llg;l, £ 2 = 5461^ [32]. Using three years of WMAP data, Spergel et al. 
obtained that the Hubble constant Hq = 73 km/s Mpc, the baryonic matter 
density Qb,o — 0.0222/i -2 , and the matter density fi mi o = 0.128/t, -2 [33], which 
are in good agreement with the observation of position of the first peak (see 
Fig. 2) but lead (assuming the ACDM model) to a value f2 m>0 = 0.24. There 
is also disagreement between H values obtained from SNIa and WMAP. We 
compute the location of the first peak as a function of Qdr,o assuming Hq = 65 
km/s Mpc (n mj o = 0.3). Now we obtain agreement with the observation of the 
location of the first peak for non-zero values of the parameter f2 dr o (Fig. 2). 
We obtain -1.05 x 10~ 5 < £l drfi < -0.5 x 10~ 5 at the 95% confidence level. 
Please note that our limit is stronger than that obtained by Ichiki et al. [12], 
which provides bounds of —7.22 x 10~ 5 < f2dr,o < 0.65 x 10~ 5 (in the case of 
the BBN) and -2.41 x 10~ 5 < fi dr , < 0.62 x 10~ 5 (in the case of the CMBR). 

All these values of fidr,o are i n agreement with the result obtained from the 
combined analysis because the 2a confidence interval for this parameter ob- 
tained from the SNIa data contains the area allowed from the CMBR. While 
the combined analysis gives the possibility that f2dr,o is equal to zero, the 
CMBR location of the first peak seems to exclude this case for h = 0.65. 



3 Conclusion 



In our paper we analysed the observational constraints on the negative {l+z) A - 
type contribution in the Friedmann equation. We discussed different proposals 
for the presence of such a dark radiation term. Although it is not possible, with 
present kinematic observations, to determine the energy densities of individual 
components which scales like radiation, we show that some stringent bounds 
on the value of this total contribution can be given. Our detailed conclusions 
are the following. 
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68.3% 



95.4% 



0.50 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 

^m,0 

Fig. 1. The 68.3% and 95.4% confidence levels (obtained from the combined analysis 
of SN+RG+SDSS+CMBR) on the plane (On,o, fi A ,o)- 



Table 1 

Results of the statistical analysis of the model with the negative radiation term 
obtained from \ 2 best fit. The upper section of the table represents the constraint 
^fc,o = (flat model). 



Sample 



^A;,0 ^m,0 ^dr,0 ^A,0 



X 



SN 

SN+RG 

SN+RG+SDSS 

SN+RG+SDSS+CMBR 



0.31 0.000 0.690 297.5 

0.32 -0.001 0.681 320.4 

0.30 0.000 0.700 322.4 

0.30 0.000 0.700 322.5 



SN 

SN+RG 
SN+RG+SDSS 



-0.28 0.43 0.000 0.850 296.0 
-0.20 0.40 -0.001 0.801 319.5 
0.06 0.29 0.000 0.650 321.1 



SN+RG+SDSS+CMBR 0.00 0.30 0.000 0.700 322.5 
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Fig. 2. The location of the first peak as a function of — f^o- Note that l\ ~ 220 for 
h = 0.73 favours f2dr,o — 0, while for h = 0.65 lead to fidr,o 7^ 0. 

Table 2 

Results of the statistical analysis of the model with the negative radiation term. 
The values of the model parameters are obtained from the marginalized likelihood 
analysis. We present the maximum likelihood values with 68.3% confidence ranges. 
The upper section of the table represents the constraint J\,0 = (flat model). 



Sample 




^m,0 


^dr,0 


^A,0 


SN 




u -' yt -0.03 


0.000_o.066 


0.66l°;° 2 


SN+RG 




S5+ - 10 


0.000_o.070 


0.65l°;° 2 


SN+RG+SDSS 




u - Ji -0.02 


0.000_o.oo7 


0.691°;°! 


SN+RG+SDSS+CMBR 




o 3n+ - 02 

u - JU -o.oi 


O.OOOOO-0.00022 


0.70l°;° 2 


SN 


-n 44+ - 24 


o.63i8:il 


O.OOO-o.iso 


o-87±8:« 


SN+RG 


U.<JO_ 2 2 


0.591°;?° 


O.OOO-o.isi 


0.82l°;i 2 


SN+RG+SDSS 


11 +°-° 8 

u - L1 -0.07 


0-29±g;g 


0.000_o.023 


61+ 05 


SN+RG+SDSS+CMBR 


00+ 02 
u - uu -0.03 


0.301°;°? 


0.00000_ .ooo26 


70+ - 02 

u - ' u -0.02 



(1) The combined analysis of SNIa data and FRIIb radio galaxies using 
baryon oscillation peaks and CMBR "shift parameter" give rise to the 
almost flat universe with Q m — 0.3. 

(2) From the above-mentioned combined analysis, we obtain an upper bound 
— fidr.o < 0.00035 at the 95% confidence level. This is a stronger limit than 
obtained previously by us from SNIa data only [21]. 
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(3) We find new stringent limits on a negative component scaling like ra- 
diation from the location of the peak in the CMBR power spectrum, 
-1.05 x KT 5 < fi dr>0 < -0.5 x 10~ 5 at the 95% confidence level. This 
bound is stronger than that obtained from the BBN and CMBR by Ichiki 
et al. [12]. 

(4) From the limit — f2dr,o < 1-05 x 10~ 5 we obtain that f2 dr ,o + ^r,o > 0. 
This implies that H 2 (z) is always greater than zero (H 2 (z) > 0) and 
the bounce does not appear which means that the big bang scenario is 
strongly favoured over the bounce scenario. 

(5) The discussed model with a small contribution of dark radiation type can 
also resolve the disagreement between values of Ho obtained from SNIa 
and WMAP. 
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